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ABSTRACT 

We report the detection of several molecular gas-phase and ice absorption features in three photometrically- 
selected young stellar object (YSO) candidates in the central 280 pc of the Milky Way. Our spectra, obtained 
with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, reveal gas-phase absorption from 
CO2 (15.0p,m), C2H2 (13.7/xm) and HCN (14.0/im). We attribute this absorption to warm, dense gas in massive 
YSOs. We also detect strong and broad 15pm CO2 ice absorption features, with a remarkable double-peaked 
structure. The prominent long-wavelength peak is due to CHsOH-rich ice grains, and is similar to those found 
in other known massive YSOs. Our IRS observations demonstrate the youth of these objects, and provide the 
first spectroscopic identification of massive YSOs in the Galactic Center. 
Subject headings: infrared: ISM — ISM: molecules — stars: formation 



1. INTRODUCTION 

The Central Molecular Zone (CMZ) is the innermost ~ 
200 pc region of the Milky Way Galaxy. It is a giant molec- 
ular cloud complex delineated by a gradient in the CO col- 
umn density and temperature. The CMZ contains ~ 10% of 
the Galaxy's molecular gas, and produces 5%-10% of its in- 
frared and Lyman conti nuum luminosities (see a review by 
iMorris & SerabyrJl 19961 and references therein). 

Evidence is mounting that conditions for star formation in 
the CMZ are significantly different from those in the Galac- 
tic disk. The gas pressure and temperature are higher in the 
CMZ than in the average disk, conditions that favor a larger 
Jeans mass for star formation and an initial mass function bi- 
ased towards more massive stars. Furthermore, the presence 
of strong magnetic fields, tidal shear, and turbulence chal- 
lenges the standard paradigm of slow gravitational collapse 
of molecular cloud cores. 

The CMZ provides several signposts of in situ star forma- 
tion, such as H2O masers, (ultra-)compact H II regions, young 
OB stars, and young supernova remnants. However, young 
stellar objects (YSOs or protostars), which are the direct trac- 
ers of current star formation, have so far eluded detection in 
the CMZ. They have be en inferred to be present based on in- 
frared photometry (e.g., 



Felli et dj2002tlSchuller et alj2006t 



lYusef-Zadeh et al]l2009l) . but spectroscopic observations are 
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required to confirm their status as a YSO. This is because 
evolved stars can look like YSOs in br oad-band photometry, i f 
they are heavily dust attenuated (e.g., Schult heis et alJl2003h . 
a problem towards the Galactic Center (GC), where Ay ~ 30. 

In this Letter, we present spectroscopic follow-up obser- 
vations of YSO can didates in the CM Z, using the Infrared 
Spectrograph (IR S: iHouck et all 120 04) onboard the Spitzer 
Space Telescope dWerner et alj|2004l) . Massive YSO candi- 
dates were photometricall y selected from the point source cat- 
alog (iRamfrez et alJl2008|). whic h was extracted from images 
of the CMZ dStolovv et aHl2006l) mad e using the Infrared Ar- 
ray Camera (IRAC; iFazio et alll2004h . This high sensitivity 
and high spatial resolution image has led to a better identifi- 
cation of YSO candidates and their follow-up spectroscopic 
observations. 

2. PHOTOMETRIC SAMPLE SELECTION 



The IRAC point source catalog (Ramirez et al. 2008) con- 
tains photometry for more than a million point sources in 
the entire CMZ (2° x 1.4° or 280 x 200 pc) in four chan- 
nels (3.6/im, 4.5/im, 5.8/j,m, and 8.0/im). Initially, we se- 
lected point sources with [3 .61 -[8.01 > 2 0, correspo nding 
to YSOs with > 2.5M Q (IWhitney et alJl200ll2004T) . We 
further confined the sample to those within \b\ < 15', result- 
ing in 1207 objects. When we had photometric measure- 
ments in at least 5 b andpasses from IRAC, 2MASS (JHK S ; 

" 2006b . and/or ISOGAL (7 pm and 15/xm; 

selected YSO candidates by com- 



Skrutskie et al 
Omont et all 120031) . 



paring the obse rved spectral energy distribution (SED) with 
Y SO models dRobitaille et al.ll2006[) using a SED fitting tool 
bv lRobitaille et all ( 120071)7 Otherwise , we ap plied additional 
color constraints from IWhitney et all (120041 [3.6] -[4.5] > 
0.5, [4.5] -[5.8] > 0.5, and [5.8] -[8.0] > 1.0) to identify 
YSO candidates. SED fitting and color selection narrowed 
down our sample to about 200 objects. 

Then, we carefully inspected IRAC three-color images to 
select objects that are distinct within the IRS slit entrances 
against the crowded stellar field and bright local background. 
Finally, a literature search was carried out for the selected ob- 
jects, and one Wolf-Rayet star and four OH/IR stars were dis- 
carded. Our final sample is composed of 107 objects, among 
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which 25 were previo usly known YSO candidates from ISO- 
GAL dFelli et al]l2002l) . 

3. IRS OBSERVATIONS AND DATA REDUCTION 

We obtained spectroscopic data for 107 YSO candidates us- 
ing the four IRS modules in May and October 2008. We ob- 
served each target in IRS staring mode with 4 exposures per 
source (2 cycles). Exposure times were 6 sec-120 sec in SH 
(short-high; short wavelength, high resolution), 6 sec-60 sec 
in LH (long-high), 6 sec-14 sec in SL (short-low), and 6 sec in 
LL (long low) modules, depending on the source's brightness, 
to achieve a signal-to-noise ratio (S/N) of at least 50 in SH and 
SL, and a minimum S/N of 10 in LH and LL. We reduced the 
IRS spectra from the basic calibrated data (BCD) products 
version S17.2.0 and S18.1.0, using the SSC software pack- 
ages IRSCLEAN (to correct for bad pixel values) and SPICE 
(to extract spectra). 

Because the GC exhibits strong, spatially variable back- 
ground, we observed multiple off-source measurements (one 
cycle, lxl mapping mode) to derive backgrounds near each 
of our YSO candidates in the four IRS modules. The on- 
source and the off-source observations were taken consecu- 
tively to minimize zodiacal light and instrumental variations. 
For the high resolution observations, we observed and ex- 
tracted four background positions (~ ±1' offsets in either 
RA. or Decl.). For the low resolution observations, we took 
spectra from two background positions at ~ ±1' away in the 
direction perpendicular to the slit, and extracted two addi- 
tional background spectra at positions along the on-source 
slit. In all of the four different IRS modules, we tried to ex- 
tract the background spectra at the same position as much as 
possible, to minimize the flux difference from different mod- 
ules. 

We made an interpolation of a plane in three dimensional 
space (positions on the IRAC map and wavelength) to obtain 
a background spectrum at the source position. We estimated 
an error in each source's background from the dispersion of 
four different background spectra, constructed from alternate 
sets of three out of the four background pointings. 

A complete analysis of spectra for all of our 107 YSO 
candidates will be presented elsewhere (D. An et al. 2009, 
in preparation). For the current analysis, we selected 
three targets (Table [U from among those showing char- 
acteristic spectral features of massive YSOs, which in- 
clude gaseous molecular absorptions from C7H7, HCN, CO? 
(e.g., Lahuis & van Dishoeck 2000; Boonmanetal. 2003; 
iKnez et aLll2009l) . a nd a solid-phase abso rption from CO2 ice 
bending mode (e.g., Gerakine s et al.lfl999l) . 

Both SSTGC 797384 and SSTGC 803187 are associated 
with a relatively weak radio c ontinuum source (SGR B2(P) 
and SGR B2(R), respectively; iMehringer etalll 19931) . They 
are on the outskirts of the Sgr B2 molecular cloud (~ 2 pc- 
4 pc from the well-studied radio source SGR B2(M)), which 
is one of the most activ e complexes of c ompact H II re- 
gions in the Galaxy (e.g., Mehringerl l 19951) . Mehringe r et al.1 
( 1993) derived zero-age main-sequence spectral types of B0 
and 06.5 for these compact H II regions, respectively, from 
the number of ionizing photons. SSTGC 524665 does not 
have radio continuum emission associat ed with it. How- 
ever, it is coincident with an H2O maser (Forst er & CaswelH 
Il989h. and is adjacent to a region of 4.5/^m excess emission 
(lYusef-Zadeh et aLll2009l). possibly tracing shocked molecu- 
lar outflows (e.g.. lSmith et al.ll2006l) . 

For SSTGC 803187, we used a non-standard extraction 



aperture in SL, because of a nearby source (« 7" south of the 
target) along the slit. We followed the prescription on the IRS 
data reduction website 11 to calibrate the flux. We trimmed 
the end of the orders to remove the noisy part of spectra, 
and spectra from different orders in high-resolution modules 
were averaged using a linear ramp. After background sub- 
traction, the SH and LH spectra were scaled down in flux to 
LL over the common wavelength interval for SSTGC 797384 
and SSTGC 803187. The SL spectra were then scaled to SH. 
For these sources, we assumed that the flux mismatch is due 
to narrower slit entrances in SH and SL. For SSTGC 524665, 
we used the SL as a basis for the scaling, because our observa- 
tions in LL and LH were contaminated by extended emission 
from a nearby (« 10" sou thwest of the target) bright source on 
the 24 /xm MIPS image dCarev et al.ll2009h I Yusef-Zadeh et alj 
2009). The background for this target is likely to be over- 
subtracted, because the target lies on a dark cloud with high 
extinction, while background spectra were taken at brighter 
spots. The potential problem of the background subtraction 
results in H2 lines (arising from the surrounding sky) appear- 
ing in absorption in SSTGC 524665. In the following initial 
analysis, we did not use LH data for all targets, but focused 
on the spectral features in other modules. 

4. ANALYSIS AND RESULTS 

Figure Q] displays background-subtracted spectra of 
SSTGC 797384, in SL (A < 11.2^m), SH (11.2/zm < 
A < 19.3/im), and LL (A > 19.3/im). The observed spec- 
trum is characterized by an extremely red SED [a = 
dlog(\F\)/dlog(X) 2], strong and deep silicate absorptions 
at 9.7 fim and 18/im, ice absorption features at 6/im, 6.85/im, 
13/im, and 15.2/im. Although the presence of forbidden lines 
indicates that these objects are likely associated with an (ultra- 
compact H II region, it could be also due to under-subtracted 
emissions from the background. 

Figure [2] shows gas-phase molecular absorptions at 
13.71/im (C 2 H 2 v s = 1- 0), 14.05 ^m (HCN v 2 = 1- 0), and 
14.97/im (C0 2 Vz = 1 - 0), detected in three YSO candi- 
dates. To derive the excitation temperature (r ex ) and col- 
umn density (N cn \) for each mo lecular species, we used model 
spectra fromlCami et alj (2006) based on HITRAN0 4 linelist 
dRothmanetalJ 120051) for C2H2 an d HCN, and those based 
on HI TEMP (lRothmanetalJll996l) for C0 2 . A second or- 
der polynomial was used to set a local continuum at 13.30/im 
< A < 14.55^m for C 2 H 2 and HCN, and 14.77/im < A < 
15.06/im for CO2. We did not include isotopes in the com- 
putation because of the limited parameter span in the model 
grids. However, ev en a relatively h igh fraction of isotopes in 
GC ( 12 C/ 13 C ps 23; IWannieri[l980l) has a negligible impact in 
the model fitting. 

We first made a fit to C2H2, and subtracted its contribution 
to the absorption near weaker HCN bands. Best-fitting model 
r ex and N co \ were found by searching the minimum x 2 of the 
fits over 100 K < T ex < 1000 K in steps of Ar ex = 100 K, and 
15 < logAU < 18 for C 2 H 2 , 16 < logAU < 18 for HCN, 
and 16 < log/V co i < 22 for CO2 with intervals of 0.1 dex. 
Solid lines in Figure|2]show our best-fitting models, and their 
r ex and N co \ are listed in Table [T] Errors in these parame- 
ters were estimated from A\ 2 , where la measurement errors 
were taken from the scatter of flux in the spectra. Systematic 
errors from background subtraction and nodding differences 
were then added in quadrature. We tested with varying cov- 

11 See http://ssc.spitzer.caltech.edu/IRS/calib 
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Wavelength (^m) 

FIG. 1.— Composite IRS spectrum of SSTGC 797384. The spectrum is from SL at A < 11.2/im, SH at 11.2/im < A < 19.3/^m, and LL at A > 19.3^tm. 
This composite spectrum is characterized by an extremely red spectral energy distribution, strong and deep silicate absorption, and several molecular gas- and 
solid-phase absorptions. 



TABLE 1 
Properties of the Sample 



Quantities 


Units 


SSTGC 524665 


SSTGC 797384 


SSTGC 803187 


R.A.(J2000.0) 


h:m:s 


17:45:39.86 


17:47:23.68 


17:47:26.29 


Decl.(J2000.0) 


d:m:s 


-29:23:23.4 


-28:23:34.6 


-28:22:1.5 


ukidss r 


mag 




18.23 ±0.06 


17.39 ±0.03 


UKIDSS H a 


mag 




14.68 ±0.01 


16.60 ±0.05 


UKIDSS 


mag 


15.71 ±0.10 


12.92 ±0.01 


14.37 ±0.02 


IRAC [3.6] b 


mag 


11.42 ±0.01 




12.22 ±0.02 


IRAC [4.5] b 


mag 


8.63 ±0.01 


9.41 ±0.01 


8.97 ±0.01 


IRAC [5.8] b 


mag 


7.08 ±0.01 


7.66 ±0.01 


7.24 ±0.01 


IRAC [8.0] b 


mag 


6.13±0.01 


5.64 ±0.01 


5.11 ±0.01 


MIPS [24] c 


mag 


1.54±0.01 


0.55 ±0.01 




Tex (C 2 H 2 gas) 


K 


300 ±150 


200 ±150 


300 ±150 


T ex (HCN gas) 


K 


400 ± 50 


100 ± 50 


100 ±50 


T ex (C0 2 gas) 


K 


200 ±50 


100 ± 50 


100 ±50 


N C oi (C 2 H 2 gas) 


10 16 cm" 2 


7.9±3.3 


1.0 ±0.3 


2.0 ±0.8 


N col (HCN gas) 


10 16 cnr 2 


15.8±4.6 


1.0 ±0.4 


2.0 ±0.9 


N col (C0 2 gas) 


10 16 cnT 2 


20.0 ±5.6 


5.0±1.5 


7.9 ±2.8 


Ncol (C0 2 solid) 


10 19 cnT 2 


0.11 ±0.01 


0.13±0.01 


0.21 ±0.01 


N co l (H 2 solid, 6/im) 


10 19 cm- 2 


< 1.7 


< 2.3 


<4.7 


Ncol (H 2 solid, 13/^m) 


10 19 cm- 2 


0.6 ±0.4 


1.3±0.4 


1.9±0.3 


Ncol (H 2 ) 


10 22 cnT 2 


2.3 ±0.2 


5.0 ±0.5 


5.8 ±0.7 


Nc 2 h 2 /Nh 2 


io- 7 


34.3 ±14.7 


2.0 ±0.6 


3.4± 1.4 


Nhcn/Nh 2 


lO" 7 


68.7 ±20.9 


2.0 ±0.8 


3.4± 1.6 


Nc0 2 ,gas/N H2 


ur 7 


87.0 ±25.5 


10.0 ±3.2 


13.6±5.1 


N C0 2 .gas/Nco 2 ,solid 
N C0 2 .solid / N H 2 0,solid 
Av (T9.7 Mm ) 




0.18 ±0.05 


0.04 ±0.01 


0.04 ±0.01 




0.18±0.12 


0.10±0.03 


0.11 ±0.02 


mag 


27 ±3 


53 ±5 


62±7 


Ay (color) d 


mag 


29.0 ±2.6 


32.0 ±2.6 


27.0±5.4 



NOTE. — Random errors are shown in the photometry. 

a A perture3 magn i tudes from UKIDSS DR2 fWarren et all l2007l) . b Photometry 
from IRamfrez et all FMM ). C Photometry from MIPSGAL (S. Carey, 2008, private 
communication). d Based on the 2MASS and IRAC color-m agnitude diagrams of GC red 
giant branch stars within 2' of the source I Schultheis et al. 2 0091) . 
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FIG. 2. — Gas-phase molecular absorptions from C2H2 1/5 = 1 — 
(13.71/xm), HCN f 2 = 1-0 (14.05/an), and C0 2 v 2 = 1 -0 (14.97/an). Best- 
fitting models are shown in solid lines. 



ering factors, but found that best-fitting case yields its value 
equal to or close to unity. 

These gaseous bandheads have been detected in absorption 
toward YSOs, tracing the warm and dense gas in the circum- 
ste liar disk and/or envelopes (e.g.. lLahuis & van Dishoeckl 
I2000t iBoonman et all 120031; iKnez et alj 120091) . They are 
sometimes detected in the photosphere and/or the circumstel- 
lar envelope of c arbon-rich asymptotic giant branch stars (e.g. 
lAoki et aTll 19991), but carbon stars have not been found in the 
GC region (e.g.. lGuglielmo et aljfl998b . 

The above estimates are based on models with a Doppler 
parameter b = 3 kms" 1 . The line width measurements of 
these molecules for several massive YSOs and that of the 
strongest H2CO absorption components ne ar SSTGC 803187 
are in the range of b = 1-7 km s" 1 (e.g. , IMehringeri ll995t 
Ivan der Tak etail 120001: IKnez et all 120091) . There are lim- 
ited model grids at b = 10 km s _1 for C2H2 and HCN, but 
Tex and N co \ were generally found within 2er from those at 
b = 3 km s" 1 . 

Figure [3] shows optical depth spectra of our sources (grey) 
at ~ 15.2/im, where the strong and wide CO2 ice absorption 
is seen. We set a local continuum over 14.5/im < A < 16.5/im 
using a 3rd order polynomial, and followed the prescrip- 
tion in Pontopp idan et al.1 (120081) to decompose the absorp- 
tion profile with five laboratory spectral components: po- 
lar (C0 2 :H 2 = 14 : 100 at 10 K; dotted line, centered at 
- 15.3/im), apolar (CO:C0 2 = 100 : 70 at 10 K; dotted line, 
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FIG. 3. — Optical depth spectra of solid-phase absorption from the CO2 ice 
bending mode. Best-fitting CO2 ice models and individual CO2 ice compo- 
nents are shown for each target: polar (dotted line, centered at ~ 15.3/im), 
apolar (dotted line, centered at ~ 15.1/xm), pure (blue shaded), diluted (black 
solid line), 15.4/xm shoulder (orange-shaded), and the sum of these absorp- 
tion components (green line). The bottom panel shows a comparison of the 
ice absorption profile between our sources (grey) and massive YSO W33A 
(blue). The optical depths for our targets were scaled in the bottom panel for 
comparison. 



centered at ~ 15.1/im), pure CO2 (15 K; blue shaded), di- 
luted C0 2 (CO:C0 2 = 100 : 4 at 10 K; black solid line), and 
15.4/im shoulder CO2 ice profile (modeled with two Gaus- 
sians in wavenumber space; orange shaded). We found a 
best-fitting set of mod els from the non- linear least squares 
fitting routine MPFIT (lMarkwardtll2009l) . Green solid line 
represents the sum of all of the ice components, and the 
CO2 ice column density in Table Q] was estimated from the 
integrated absorption, adoptin g the integrated line strength 
A = 1.1 x 10" 17 cm molecule" 1 fcerakines et alj|1995l) . 

Unlike the CO2 abs orption profiles obs erved in quiescent 
molecular clouds (e.g. lWhittet et alj|2009l) . the 15.2//mband 
in Figure [3] shows a remarkable double-peaked profile. Dou- 
ble peaked profiles are commonly observed towar d YSOs 
(e.g jGerakines et aTll 1 9991 ; IPo ntoppidan et al. 2008), and are 
ascribed to pure CO2 ices resulting from crystallization of 
heated H20-rich ices. However, the double peaks toward 
the GC candidate YSOs are centered at longer wavelengths 
(15.15/imand 15.4/imvs. 15.10/imand 15.25/im), and result 
from CO-rich (15.15/impeak) and CHbOH-rich ices (15.4/im 
peak; see Fig. [3]). 

The strength of the 15.4/im peak is similar to that of the 
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Best-fitting parameters are listed in Table Q] We ob- 
tained a total hydrogen column density from the optical 
depth of the 9.7/im silicate absorption, assuming Ay/rgj = 9 
dRoche& Aitkenl 19851) and Ah/Ay w 1.87 x 10 21 cm^mag" 1 
Bohli n et alJI 19781) at Ry = 3.1. The H2 column density was 



a 
a 

o 

a 

o 



8 910 20 
Wavelength (yu.m) 

FIG. 4.— Fit to the H 2 ice and silicate absorption for SSTGC 797384. 
Top: SL and LL data (grey), with a best-fitting pseudo-continuum (black 
line). Bottom: Decomposition of optical depth spectra (grey) with the silicate 
(red) and the laboratory FLO ice profiles (blue). Black line represents a sum 
of these two components. 



well-studied embedded massive YSO W33A (Gera kines etafl 
1999, bottom panel in Fig. [3J. It is ascribed to a Lewis 
acid-base interaction of CO2 (the Lewis acid) with CH3OH 
dDartois et alj[l99 9a). Other species could be acting as a base 
as well, but CH3OH is preferred due to its high abundance 
toward W33A: 5%-22% relative to solid H9O dDartois etalj 
1999b). Two other YSOs (AFGL 7009S, AFGL 2136) show a 
prominent 15.4/im peak, and i ndeed these sources have high 
CH3O H abundances as well dDartois et alJll999btlGibb et alj 
120041) . This suggests that the GC candidate YSOs have high 
solid CH3OH abundances as well. 12 Although the origin of 
the large quantities of CH3OH in th e previously studied m as- 
sive YSOs is not fully understood (Da rtois et al.lll999a) . so 
far all lines of sight with high solid CH3OH abundances are 
associated with star formation, strengthening the idea that the 
sources studied in this paper are indeed YSOs. 

To derive abundances of these molecular absorptions with 
respect to the hydrogen a nd solid H2O column densities, we 
followed the procedure in Boogert et al. (2008) to fit the H2O 
ice and silicate absorption profiles to SL and LL spectra. 
Figure shows an example for SSTGC 797384. We used 
the silicate absor ption profiles in the line of sight to the GC 
(GCS 3 spectrum; Kemp eret alj2004l) plus a la boratory spec- 
trum of pure amorphous H2O ice at T = 10 K (Hudgi ns et alj 
119931) . We simultaneously fit a second-order polynomial for a 
pseudo-continuum (i.e., including corrections for the continu- 
ous extinction), the silicate profile, and H2O ice absorption to 
the 5/im < A < 32/im spectrum. We masked absorption fea- 
tures at 6/xm, 7/xm, and 15/im, and all unresolved emission 
lines, before performing a non-linear least squares fit. 

12 This needs to be verified by indepen dent L-band spectro scopy of the 
3.53/mi C-H stretch mode of CH3OH fe.g JDartois et aU1999rj) . 



then approximated by 7Vh 2 = Nu/2. The ice column density 
for the 13/.tm librational H2O absorption was estimated from 
the integrated absorption of the best-fitting H2O model. The 
H2O ice column density from the 6/im bending mode, fit 
separately after fixing the continuum and extinction to pre- 
viously found values, is an upper limit because the 6/im ab- 
sorption is not solely due to H2O ice. We adopted the in- 
tegrated line strengths A = 1.2 x 10~ 17 cm molecule" 1 for the 
bending mod e andA = 3.1 x 10~ 17 cm molecule" 1 for the libra- 
tional mode (Gerakines et al. 1995). Errors in these parame- 
ters (Table [U are formal estimates made by varying the range 
of wavelengths that we used for the 9.7 /jm silicate fitting, or 
by taking a few different ways of setting the continuum. 

The gas-phase molecular abundances relative to H2 are 
listed in Table Q] Our derived abundances of ~ 10~ 7 -10~ 6 
for C2H2 a nd HCN are comparable to those found for mas - 
sive YSOs dLahuis & van Dishoeckll2000l: iKnez et alj|2009b . 
although abundances for SSTGC 524665 have large errors. 
Intervening molecular clouds in the line of sight to the GC 
are less likely the main cause of these abs orptions, because 
the av erage HCN abundance of 2.5 x 10~ 8 dGreaves & Nvmaril 
1996) towards Sgr B2(M) is an order of magnitude lower 
than our measurements. Our gas-phase CO2 abundances are 
an order of ma gnitude larger t han those found towards mas- 
sive YSOs in iBoonman et alj d2003l) . but our gas to solid 
abundance ratios for CO2 are consistent with their estimates 
(10 _1 -10~ 2 ). Our abundance of CO2 ice relative to H2O ice 
is within the range (0. 10-0.23) found towards massive YSOs 
dGerakines et alj|1999l) . 

Finally, Table Q] lists our estimate s on Ay from the 9.1 jjtm 
silicate absorption and those from Schulth els et al.1 (2009), 
based on the 2MASS and IRAC color-magnitude diagrams 
of GC red giant branch stars within 2' of the source. Both 
SSTGC 797384 and SSTGC 803187 have higher Ay val- 
ues than the average for field stars, implying that a signif- 
icant fraction of the attenuation is intrinsic to the source. 
SSTGC 524665 has a lower Ay, comparable to the aver- 
age value for surrounding field stars. We also note that 
SSTGC 524665 is located at b » -0.2°, so it is possible that 
it is in front of the GC. If we assume a distance of 8 kpc for 
all three sources, and adopt the extinction of surrounding field 
stars as the foreground extinction to each source, then we de- 
rive stellar masses of 12±3M , 14±3M , 17±6M for 
SSTGC 524665, SSTGC 797384, and SST GC 803187, re- 
spect i vely, b y using a grid of YSO models (Robitaille et al. 
2006, 2007). More detailed discussion of the model fitting 
will presented in a future paper. 

To summarize, we presented the evidence from IRS spec- 
tra for the first spectroscopic identification of massive YSOs 
in the GC. In our next paper (D. An, 2009, in preparation), 
we will present the results for all 107 YSO candidates, to- 
gether with additional data from millimeter to radio observa- 
tions, and use them to better understand the nature of these 
embedded sources. 
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